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Abstract 


In  this  paper  we  develop  local  and  global  estimates  for  the  solution  of 
convection-diffusion  problems.  We  then  study  the  convergence  properties  of 
a  Time  Marching  Algorithm  solving  Advection-Diffusion  problems  on  two 
domains  using  incompatible  discretizations.  This  study  is  based  on  a  De- 
Giorgi-Nash  maximum  principle. 
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1  Introduction 


Domain  decomposition  methods  have  recently  become  an  eiRcient  strategy 
for  solving  large  scale  problems  on  parallel  computers  ([1],  [2],  [3],  [4],  [5], 
[6]).  Nevertheless,  they  can  also  he  used  in  order  to  couple  diiferent  models 
[11],  [18].  In  this  paper  we  will  examine  a  domain  decomposition  strategy 
which  can  be  apphed  to  such  situations. 

This  approach  was  introduced  in  order  to  solve  several  difficulties  that 
occur  in  fluid  mechanics.  In  particular,  our  aim  is  to  introduce  several  sub- 
domains  in  order  to  do  one  of  the  following  : 

•  Solve  dilferent  problems  on  each  subdomain. 

•  Use  different  kinds  of  approximation  methods  on  each  subdomain  [7]. 

•  Use  “local  refinement  techniques”  or  “mesh  adaptive  techniques”,  lo¬ 
cally,  per  subdomain  ([10]). 

The  subdomains  fuUy  overlap  and  the  coupling  is  achieved  through  “fric¬ 
tion”  forces  acting  on  the  internal  boundary  of  the  domain,  these  friction 
forces  being  updated  by  an  explicit  time  marching  algorithm. 

Several  versions  of  this  methodology  have  been  studied  in  [15].  In  [15]  the 
empha.siz  wan  on  the  implicit  time  discretization  version  of  this  algorithm, 
we  focus  in  this  paper  on  the  explicit  version  of  this  methodology.  The 
theoretical  study  of  our  method  will  be  done  on  an  Advection-Diffusion 
problem,  which  will  serve  as  our  model  problem.  The  analysis  will  be  made 
at  the  continuous  level,  independently  of  any  (space)  discretization  strategy, 
which  means  that  the  derived  results  will  be  mesh  independent. 

In  the  next  section  we  develop  a  maximum  principle  for  general  second- 
order  elliptic  problem  based  on  the  De-Giorgi-Nash  theory.  In  section  2  and 
3,  we  develop  estimates  for  the  solution  of  the  convection-diffusion  problems 
respectively  with  Direchlet-Neumann  and  Direchlet  boundary  conditions. 
These  results  are  based  on  the  maximum  principle  of  section  2.  We  then 
apply  these  tools  to  the  analysis  of  an  explicit  time  marching  algorithm.  We 
also  study  a  fixed  point  method  for  the  implicit  time  marching  algorithm  of 
[15].  Practical  applications  of  the  time  marching  algorithm  to  real  life  CFD 
problems  can  be  found  in  [14],  [19],  [20],  and  [21]. 
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2  Local  estimates 

In  this  section  we  shah  establish  a  maximum  principle  for  an  arbitrary  elliptic 
operator  of  second  order.  These  tools  are  central  to  the  development  of 
our  theory  in  order  to  derive  the  convergence  analysis  of  the  explicit  time 
marching  algorithm  described  in  section  (5.1). 

Let  L  be  an  operator  written  under  the  form 

Lu  =  a'^{x)DijU  +  V{x)DiU  +  c(x)«, 

for  any  u  in  with  0,  a  bounded  domain  of  IR"^.  The  coefficients 

and  c,i,j  =  l,...,n  are  defined  on  Cl.  As  usual,  the  repeated  indices 
indicate  a  sximmation  from  1  to  n. 

We  suppose  that  the  operator  L  is  strictly  elliptic  in  in  the  sense  that 
the  matrix  A  of  coefficients  [a^^]  is  strictly  positive  everywhere  in  Cl.  Let  A 
and  A  denote  respectively  the  smahest  and  the  largest  eigenvalue  of  A.  Let 
V  denote  the  determinant  of  the  matrix  A  and  V*  =  We  have 

0  <  A  <  D’  <  A. 

We  suppose  in  addition  that  the  coefficients  0*^,6*  and  c  are  bounded  in  Cl, 
and  that  there  exists  two  positive  real  numbers  7  and  6  such  that  : 

<  7»  (L  is  uniformly  elliptic)  (1) 

<  s.  (2) 

Now,  we  are  in  a  position  to  state  the  principal  result  of  this  section, 
proved  in  annex. 

Theorem  2.1  Let  u  e  W^’^(Cl)  and  suppose  that  Lu  >  f  with  f  €  L^{Cl) 
and  c  <  0.  Then  for  all  spheres  B  =  B2R{y)  of  center  y  and  radius  2R 
included  in  Cl  and  for  all  p>  Q,  we  have  : 

S  fll/ll”.®}-  (3) 

where  the  constant  Cr  depends  on  {n,'f,6R?,p),  but  is  independent  of  c. 
Above  u"*"  =  max{u,  0). 
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Remark  2.1  The  statement  of  the  same  theorem  can  be  found  in  [12],  un¬ 
der  the  assumption 


|c|/A<i.  (4) 

So^  there  the  constant  Cr  depends  indirectly  on  c  through  6.  That  is 
exactly  what  we  want  to  avoid,  since  we  would  like  this  constant  to  be  inde¬ 
pendent  of  c  (see  section  5.1). 

3  First  fundamental  estimate 

Let  Qioc  Le  a  connected  domain  of  IR”  with  ilioc  C  (Figure  1).  The 
boundaries  of  the  two  subdomains  are  defined  as  foflows: 

Tj  =  do,  n  d^lioc,  ( internal  boundary) 

Fj  =  dQ,ioc  n  (  interface) 

Too  =  (farfield  boundary) 

We  denote  by  n  the  external  unit  normal  vector  to  dO,  or  dQ,ioc- 

Let  y  be  a  given  velocity  field  of  an  inviscid  incompressible  flow  such  that: 

(  divV  =  0  in 

{  (5) 

V.n  =  0  on  r^. 

We  shall  derive  an  estimate  for  the  solution  of  the  following  Direchlet- 

—i/Av  +  V  •  Vu  +  —V  in  ft,  (6) 

r 

0  on  Too?  (7) 

^  on  (8) 

where  the  function  g  is  given  in  and  the  coefficient  r  is  strictly 

positive,  and  u  is  the  diffusion  coefficient.  Let  W  be  the  sub-space  of 
defined  by 


Neumann  problem: 


Cv  = 

V  = 

dv  __ 
dn 


W  =  {w^  u;  =  0  on  Too} 


(9) 
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We  then  define  the  following  bilinear  forms  on  W 


a(v^w)=  1  uVvVw+  / 

•I  Cl  J  m 

(10) 

(v,  w)  =  f  vw. 

(11) 

The  first  basic  problem  associated  to  (6)-(8),  can  be  written  as  follows:  Find 
V  eW  satisfying  : 

a{v,w)  +  {\lT){v,w)=  I  gwdT,  Vw  £  W,  (12) 

Moreover,  we  assume  that  the  coefficients  1/  and  r  satisfy  the  following 
relation: 

I'T  <  1-  (13) 

This  hypothesis  is  not  necessary  but  simplifies  the  proofs  to  come.  More¬ 
over,  it  is  not  restrictive,  since  we  would  like  the  convergence  for  small  r 
(see  section  5.1). 

Let  d  denotes  the  overlapping  distance  as  described  in  the  Figure  2.  Let 
then  /5  be  a  real  number  such  that 

0  <  /?  <  Z^/u/d, 

and  set 

k  = 

The  first  basic  result  states  the  global  estimate  of  the  solution  of  the 
first  basic  problem  (12)  in  terms  of  the  boundary  data  g. 

Lemma  3.1  There  exists  a  constant  Co  such  that  we  have: 

Ill’ll!, <  (Co/l^)||5||-l/2,rt,  (14) 

Proof  of  lemma  3.1 

By  using  the  relation  (5)  we  have  the  following  equality: 


J  udiv(Pu) 


0,  V  u  €  W. 
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(15) 


Choosing  ty  =  in  (12),  we  then  obtain 

f  {v\Vv\^  +  {\It)v^}  =  f  gv. 

Jq  JTt 

From  this  equality  we  deduce  the  following  estimate: 

^ll^lli.n  ^  lbll-i/2,rjl^lli/2,rfc- 

The  application  of  the  trace  theorem  yields  the  estimate  (14),  which  implies 
in  particular 


lb||o,n  <  (co/i')||5||-i/2A.  (16) 

■ 

Let  flj  be  the  subdomain  of  width  |  with  external  boundary  Fj  as  described 
in  the  Figure  2.  Let  Ky  =  Bd{y)  be  the  sphere  of  center  y  and  radius 
There  exist  yi,..-,yi  belonging  to  fl,-  such  that 

^2i  =  U3,eJ2,5a(y)  C  U^j^^Ky.. 

6 

We  define  then  K  by  setting 

K  = 

The  next  lemma  states  the  local  estimate  of  the  solution  v  of  the  first  basic 
problem  (12). 

Lemma  3.2  There  exists  a  constant  ci  such  that: 

Halloo, <  ci||u||o,n.  (17) 

where  ci  is  a  constant  depending  only  on  u,‘y,8d?  and  (3/2(i)^/^. 

Proof  of  lemma  3.2 
The  operator 

L  =  -C 

satisfies  the  assumptions  of  the  theorem  2.1,  with  c  =  — l/r  and  /  =  0. 
Applying  then  this  theorem  with  p  =  2,  y  €  fij-  we  obtain 

||'*^||oo,K’y  <  ^lll^llo,B2d/3(y)' 
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Therefore 


(18) 


lbllco,ji:3,  <  ci||t7||o,n, 

where  ci  is  a  constant  depending  only  on  v,‘y,6d?  and  (3/2^)"/^. 

Applying  the  relation  (18)  to  each  Ky.  we  obtain: 

||^^lioo,/f  <  sup  Cij||u||o,fi. 

Setting  Cl  =  sup  cij,  we  finally  have 

lblloo,7r  <  ci||u||o,fi.  (19) 

And  the  lemma  is  proved. 

■ 

We  shall  now  establish  other  local  estimates  for  the  solution  v  of  the  first 
second  basic  problem.  For  any  Mi  in  fl,-,  we  introduce  (see  Figure  2): 

•  Bi  —  the  ball  centered  on  Mi  of  radius  d/6, 

•  Vi  =  exp[k(r^  -  dV36)]||u(|oo,9B,. 

We  then  have: 

Lemma  3.3  The  solution  v  of  the  first  basic  problem  satisfies: 

|u(Mi)|  <  exp(-A:dV36)||u(|oo,9Si,  VM,-  €  fi,-.  (20) 

Proof  of  lemma  3.3 

The  operator  C  applied  to  u,-,  can  be  written  in  polar  coordinates  (with 
r  =  MiM)  : 

Lvi  =  4{-k^i/r^  -  ki/  +  ^V.Crr  +  — )u,-. 

2  4r 

Therefore 

Lvi  >  4{-k^nr^  -  ^\V.er\r  +  (^  -  kv))vi.  (21) 

We  set  then  : 

<p{r,  k)  =  a(fc)r2  +  b{k)r  +  c{k),  (22) 

with 
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a{k)  =  —k^u 
m  =  -^\V.er\ 

c(k)  =  -p —  ku. 

At 

We  seek  to  satisfy  tke  followiag  relation  : 

0  <  inf  ip(r,  k)  for  0  <  r  < 

6 

As  fp{T,k)  decrecises  on  JR'*',  this  will  be  satisfied  iff 

(p{d/&)  >  0, 

i.e.  iff 
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In  addition,  by  construction 


Vi  >  V  on  dBi. 

Consequently,  by  using  the  maximum  principle  we  obtain  the  following  re¬ 
lation  : 

V  <  Vi  in  Bi . 

In  particular 

v{Mi)  <  exp{-kd? lZQ)\\v\\oo,dBi- 

We  do  tlie  same  for  —v^  and  finally  we  have 

\v(^Mi)\  <  exp(-“A?d  /36)||t;||oo,9j5^5  (24) 


Let  Qoo  be  defined  by  f2oo  =  ft  \  ft/oc-  The  next  result  establishes  an 
estimate  of  the  solution  v  of  the  first  basic  problem  on  the  domain  fto^. 

Lemma  3.4  There  exists  a  constant  C2  such  that: 

l|v||i,noo  <  M\co,ni\jc2/d  (l  -f  .  (25) 

Proof  of  lemma  3.4 
Let  f  6  be  such  that  : 


On  the  other  hand,  we  have  : 

/  (iiv(y^;)e2^  =  /  div(ye^nV2)  -  /  V.V^^v\  (28) 

Ja  JQ  JQ 

Using  the  relations  (27)  and  (28),  (26)  becomes 
0  =  jJ^(i/|V(ex;)|2  +  div(nV/2)  +  eV/r)-^(i/n2|Ve|2  +  t;2^U.V0 
=  ^  Kiv(^n)|2 + -  u)ev^  -  + v^v.vo 

=  /  i/(|Vt;p  +  |t;p)  +  [  i'(|V(^t;)|^  +  l^vp)  +  f  (1/r  — 

-  f 

J  j 

Hence,  we  obtain  : 

+  /  K|V(^n)p  +  l^np)  +  /  (1/r  -  r/)fV  = 

«/f2t  Vi* 

I  (i/t;2|vep  +  t;2eu-v0. 

jQi 

Tbe  relation  (13)  then  yields 


^U^WiQooUQi  < 

/  (^n^ivep  +  n^^y.VO 

«/f2  j 

(29) 

< 

/("ivtf +  fv.vf) 

V  S*t 

(30) 

< 

< 

(31) 

If  we  take  ^  such  that 


IlflloA  <  1, 


l^llA  “ 

where  C2  is  a  constant,  (31)  then  becomes 

||«||i,noo  <  h\\oo,ai\fc^  ,  (32) 

which  is  the  conclusion  of  our  lemma. 

■ 

Now  we  are  in  a  position  to  state  the  main  result  of  this  section. 

Theorem  3.1  Let  v  he  the  solution  of  the  first  basic  problem  (12).  If  r  is 
sufficiently  small,  we  have 


Hi/2, r,  <  + 

(l/i/)exp(-A;dV36)||5||_i/2,rfc, 

where  Ci  and  C2  are  constants,  with  C\  depending  only  on  d,  z/,7  and  S,  but 
not  on  T. 

Proof  of  theorem  3.1 

The  proof  of  this  theorem  is  based  on  the  above  lemmas.  Since  dBi  C  K, 
We  have 


Iblloo.as,-  <  ll^^lloo, iiT, 

The  lemma  3.2  then  implies 

Using  the  lemma  3.3  and  the  above  estimate  we  obtain: 

\v{Mi)\  <  ea;p(-fc</V36)ci||u||o,ji,  VM,-  €  ft,-. 
Consequently  we  have 

ll^'lloo.Qi  <  exp(-fe<fV36)ci||t;||o,fi. 


(33) 

(34) 


(35) 
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Applying  the  lemma  3.1  we  obtain: 

IbllooA-  ^  ^ca:K-^<^V36)||i7||-i/2A-  (36) 

The  application  of  the  lemma  3.4  then  yields: 

||^||i,noo  <  cociy^  (^1  +  (37) 

(l/i/)exp(-MV36)||p||_i/2A. 

To  conclude  we  use  the  trace  theorem  which  yields 

ll^lli/2,ri  <  C3||u||i,neo- 
Consequently,  we  have  the  final  estimate: 

||^’lli/2A  <  CoCiC3y^(H- ® 

(l/i')exp(-A:dV36)||sr||_i/2,r4» 

which  corresponds  to  our  theorem  with  Ci  =  C0C1C3  and  C2  —  C2. 


4  Second  fundamental  estimate 

In  this  section  we  shall  derive  an  estimate  of  the  solution  of  the  following 
Direchlet  problem: 


JOv  = 

—i/Av  -h  V  ■  Vv  -f-  —V  in  fl/oc, 

T 

(38) 

V  = 

h,  on  Tj, 

(39) 

V  = 

0,  on  Tj,, 

(40) 

where  the  function  h  is  given  in  the  coefficient  r  is  strictly  positive, 

and  V  is  the  diffusion  coefficient.  The  velocity  field  V  is  given  by  the  relation 
(5).  Let  W  be  the  sub-space  of  defined  by 
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Figure  1;  Description  of  the  Domain  Slioc  and  of  the  splitting  used  in  the 
majoration  of  the  local  solution. 
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W  =  {w  e  5^^(Q;oc)|  =  0  on  Ffc}. 

We  then  define  the  following  bilinear  forms  in  W : 


a{v,w)  =  u  /  Vv,Vw  + 

1  div(y^;)ty, 

(41) 

J^loc 

(v,w)=  [ 

j^loc 

VW, 

(42) 

^^loc 

with  V  and  w  in  W.  The  second  basic  problem  associated  to  (38)-(40) 
corresponds  to  the  following  problem; 

Find  V  such  that 

a{v,w)  +  {ljT){v,w)  =  I  v-^w,  Vn;  €  W,  (43) 

JTi  0’^ 

nir,-  =  h,  (44) 

where  h  is  given  in  We  first  have  the  following  lemma  ; 

Lemma  4,1  For  r  sufficiently  small,  we  have 

a{w,w)  +  {1/t){w,w)  >  (i'/2)||n;||f^0io.’  ^ 


Proof  of  lemma  4.1: 

Under  the  hypothesis  1/r  >  u/2  +  (l/2v)\\V\^,  and  nsing  the  Cauchy- 
Schwarz  inequality,  we  obtain  : 

a(u,u)  + (1/t)(u, u)  =  f  i/Vu.Vu+  j  V.Vvv  +  (I/t)  j 

>  +  {l/r)\\v\\l,  -  |lU||oo||Vu||o,2||u||o,2 

>  i^||Vu||2,2+(l/r)||ul|^,2-(»'/2)||Vu||o% 

>  W2)(||Vu||2,2  +  ||u||2,2)- 
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We  will  also  make  the  simplifying  assumption  (13).  We  first  establish  a 
global  estimate  for  the  solution  of  the  second  basic  problem. 

Lemma  4.2  The  solution  v  of  the  second  basic  problem  (43)-(44)  satisfies: 

IhlliAoc  2(1  +  ^  _±Hoo  j  ||h||i/2,ri  (45) 

Proof  of  lemma  4.2: 

Choosing  w  =  v  in  (43)  we  obtain  : 

u  f  |Vup+  /  (div(Vv)v  +  (l/T)v^)=  f  up-h.  (46) 

The  lemma  4.1  then  yields 

ll^lliAoc  -  2||^V^^II-i/2,rill^lli/2,r<.  (47) 

We  shall  now  establish  an  estimate  of  ||5v/5n||_i/2,r,-.  Combining  (43)  and 
(5)  we  obtain: 

f  dv  r  ,  1 

~  Jo  (^^^^  4"  (l/t')l^.Vnw  4 - vw). 

Therefore,  for  any  w  in  W,  we  have 

(livens, +  (l/r^)||«.||3,a,„)'/^ 

<  [1  +  (l±|aiy ]I/!||„||,_„,^(1 + 
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The  trace  theorem  then  yields 


lia»/s»ll-i/w  <  (1  +  A  +  ||.|b,n,„- 

(48) 

Combining  now  the  relations  (47)  and  (48)  we  have 

W|iA„  <  2(1  +  ||l.||./2,r,  (49) 

and  hence  in  particular 

ll»lloA«  <  2(1  +  ---")*'*  ( 1  +  '  ll'>lli/2,r, .  (50) 


Let  Ky  —  be  the  sphere  centered  on  y  and  of  radius  d/4,  with  y 

belonging  to  ty  (see  Figure  3).  By  construction,  Ty  is  the  center  surface  of 
Q>ioc  and  Qi  is  the  subdomain  of  width  |  centered  on  Fy  , 

We  have  the  following  lemma: 

Lemma  4.3  There  exists  a  constant  ci  such  that: 

||^^||oo,ii:  <  ci||t>||oAoc-  (51) 

Proof  of  lemma  4.3: 

Following  the  same  argument  as  in  the  proof  of  the  lemma  3.2  we  obtain: 

lkllco,iirj,  <  ci||n||oAoc>  (52) 

where  ci  is  a  constant  depending  only  on  d,  z/,  7  and  6,  On  the  other  hand 
there  exist  ^i, . . . ,  in  Qi  such  that 

l 

^2i  =  U  ^i(y)  c  U  =  K. 

j/efii  ®  3=1 

By  applying  the  relation  (52)  to  each  Kyj ,  we  obtain 
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Figure  2:  Description  of  the  local  domain  Slioc  and  of  the  splitting  used  in 
the  majoration  of  the  global  solution. 


<  sup  ci,i|t?||oAoc  =  ci||v||oAoc-  (53) 


Next  we  shall  establish  another  local  estimate  for  the  solution  of  the  second 
basic  problem  (43)-(44).  For  any  M,  €  fi,-,  we  introduce  (see  Figure  3): 

•  a  ball  Bi  centered  on  M,-  and  of  radius  d/6, 

•  the  function  u,-  =  explkij^  -  d2/36)]||n||oo,aSi. 

We  then  have: 

Lemma  4.4  The  solution  v  of  the  second  basic  problem  (43)-(44)  satisfies: 

|t7(Mi)|  <  exp[-kd? /Z%]\\v\\^^3Bi.  (54) 

Proof  of  lemma  (4.4): 

By  construction  of  k  (see  the  previous  section),  <p{r,  k)  is  positive  for  all 
r  e  [0,d/6].  Then  by  following  the  same  argument  as  in  the  proof  of  the 
lemma  3.3  we  obtain  the  inequality  (54). 
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Let  fib  be  the  subdomain  of  flioc  described  in  the  Figure  3.  The  global 
estimate  of  the  solution  of  the  second  basic  problem,  is  obtained  in  the  next 
lemma. 


Lemma  4.5  The  solution  v  of  the  second  basic  problem  (4S)-(44)  satisfies: 

<  lb||ooA\/^  (l  +  •  (55) 


Proof  of  lemma  (4*5): 

Consider  ^  such  that: 

f  ^  =  1  in 

1  supp^  C  1^6  U 

Choosing  w  =  in  (43)  we  obtain  : 


(56) 


I  {—i/Av  +  div(Fu)  +  =  0.  (57) 

f^loc 

Similarly  to  the  proof  of  the  lemma  3.4  we  obtain: 

<  /  (i/t;2|V^|2  +  t;2|y.VO. 

Choosing  ^  such  that 

ll^lloA  ^  1 

and 

Iflta  =  <^2/'*. 

we  finally  obtain  as  in  the  proof  of  the  lemma  3.4  the  inequality  (55). 

■ 

Finally,  the  main  result  of  this  section  is  presented  in  the  following  the¬ 
orem: 

Theorem  4.1  For  r  sufficiently  small,  the  solution  v  of  the  problem  (4S)- 
(44)  satisfies: 
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(1  +  l/r2)ea:p(-fedV36)||/i||i/2,rp  (58) 

where  C\  and  C2  are  constants  with  Ci  depending  only  on  d,v,u  and  6. 
Proof  of  theorem  4.1: 

Since  dBi  C  K  hy  construction,  the  lemmas  4.3  and  4.4  imply: 


ll^llooA-  e®i>(-^<^^/36)ci|lt;||o,j2,„^. 
Furthermore  by  using  the  lemma  4.2  it  follows: 


|'^||oo,f2i  ^2  1  + 


l±MloV'“ 


c,(l  +  l/r=)'/=e*p(-WV36)||4||./2,r,. 
By  using  the  lemma  4.5  we  then  obtain: 


|'^’||i,ni,uni  S  2  1  + 


i  +  l|i'll2oV'' 


1/2 


ciy/^  (1  + 

(1  +  l/T2)i/2exp(-A:dV36)||h||i/2,r,. 

Before  concluding  we  shall  establish  an  estimate  of  the  term 


||du/an||_i/2,r,. 


Choosing  w  such  that: 


w  e  n'^(Q,ioc),  with  u;  =  0  on  d^b  n  d^li, 
and  using  (43)  we  obtain: 


(59) 


(60) 


(61) 
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—i/Av  +  divCl^i;)  +  v/t)w  =  0. 


Applying  the  Green’s  formnla  and  using  (5),  we  obtain: 


f  d'O  f  1 

/  .—'w—  I  (VvVw  +  (l/i/)V.Vvw - vw). 

M  On  Juk 


Similarly  to  the  proof  of  the  lemma  4.2  we  obtain  the  following  inequality: 


l|a»/9»ll-i/2,r.  <  (1  +  (62) 

The  completion  of  the  proof  of  the  theorem  results  from  the  combination  of 
the  relation  (61)  with  (62). 


5  Convergence  analysis  of  the  explicit  time  march¬ 
ing  algorithm 

Consider  the  following  elliptic  problem: 

—  +  V.V4>  —  uA4>  =  0  in 

T 

<  <t>  =  4>oo  on  Too,  (63) 

^  =  0  on  r>, 

that  we  would  like  to  solve  by  the  fundamental  algorithm  of  [15].  This 
algorithm  can  be  written  in  this  case  as 

•  set  =  <f>oi  and  4>°  =  4>o. 

•  then,  for  n  >  Q,4>'ioc  <i>^  being  known, 
solve 

^  =  0  in  fi/oc, 

T 

'  Ct'  =  r  on  Ti,  (64) 
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I  = 

*  T 

^"+1  = 


0  in  fi, 

(f>oo  on  r  c 


(65) 


dn 


=  u 


dn 


on  Tj. 


We  shall  show  in  this  section  that  this  algorithm  converges,  and  the 
converged  solution  corresponds  to  the  solution  of  the  initial  problem  (63). 
More  precisely  we  have  the  following  theorem. 


Theorem  5.1  Forr  sufficiently  small,  (j)  being  the  solution  of  the  stationary 
problem  (63),  we  have  : 

i)  converges  to  ^  in 

ii)  converges  to  <f>  in 

iii)  converges  to  (f  in 

iv)  converges  to  4>  in  H^filioc). 


Proof  of  theorem  5.1: 

By  the  transformation  -*•  with  <f>  the  solution  of  the  stationary 

problem,  this  problem  can  be  reduced  to  the  case  (^oo  =  0.  Multiplying  the 
equation  in  (64)  by  w  6  W,  integrating  by  parts,  we  obtain: 

f  r  r 

L  L  +  ^  f  V<^7+'Vu; 

=  " /r,  ^ 

We  now  apply  the  theorem  4.1  and  we  obtain 

I -  /  -I 

Il“^-||-i/2,r6  <  ^1  +  —(1  +  ||T||^)J  (67) 
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(i + l\\v\u4^^'''‘ 

(l+l/r=)exp(-MV36)||,^“||i/,,r,. 

On  tlie  other  hand,  multiplying  the  equation  in  (65)  hy  w  and  inte¬ 
grating  by  parts  we  obtain  the  equality 

/  — - w+  f  +  u  f  —  u  f  w  (68) 

Ju  T  Jq  Ju  JTk  on 

with  w  €  H^{Q,)  and  u;  =  0  on  Too.  Applying  the  theorem  3.1  to  this 
problem  yields: 


exp(-A:dV36)||5CtV^»ll-i/2A* 

Combinig  (67)  and  (69),  we  then  have: 


(69) 


p^r.tv^^ii-i/2A 


<  c[y/^d(l  +  -^{l+\\V\\l)) 

(i+^imi=oV^)'^' 

(1  +  1/T^)exp  II^Cc/^^ll-i/2,r», 


with  k  = 


/? 

U^/T' 


Therefore  for  r  sufficiently  small,  the  coefficient  of  reduction 


will  be  domiaated  by  the  exponential  term  and  wiU  then  be  strictly  less 
than  1,  implying  the  linear  convergence  to  zero  of 


||a#r+V9»ll-i/2,r.- 
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Tlus  corresponds  exactly  to  the  statement  (i).  This  statement  combined 
with  (69)  leads  to  the  convergence  of  to  0  in  F^/2(r,).  Applying  (14) 
with  g  =  we  have  in  addition 

and  therefore  ||i,n  converges  to  zero  at  the  speed  of 
Applying  now  (45)  with  h  =  <^”,  we  also  have 

II^St'lliA^  <  2(1  + 1/’-")''*  (i  +  ^(1  +  l|i"l£,))  mWnxr 

And  then  also  converges  to  zero  at  the  speed  of  ||<^”||i/2,ri- 


5.1  Convergence  of  a  fixed  point  method  for  the  implicit 
time  marching  algorithm 

The  implicit  time  marching  algorithm  of  [15]  couples  the  global  and  the  local 
problem.  To  uncouple  them,  it  is  advisable  to  use  the  fixed  point  algorithm 
below  : 

•  set  =  V’oZ  and  (f>°  =  V’o, 

•  then  for  >  0,  being  known, 
solve 


(  •K*" 


loc,k-\-l 


<!>] 


loc 


+  <iiv(u<Af +fc+i )  -  = 


At 


jn+1 

^loc,k+l 

= 

on 

r.-, 

®Zoc,fc+l 

=  0  on 

r6, 

(70) 

Ktl  -  r 

At 


+  div(t;«^^+J)  -  = 

^k+i  - 

vd(^lfjdn  = 


0  in  ft, 

^oo  on  Tcoj 

tn+l 


on  Tj. 

(71) 
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We  will  study  now  the  algorithm  (70)-(71).  By  setting 


‘^loc,k,q  —  ^/oc,g+l’ 

(73) 

we  have  that  'ipioc,k,g  and  verify  the  following  equations  : 


f  '^loc,k,ql^t  +  div{v'tploc,k,g)  -  vA‘tj}ioc,k,q 


0  in  ^loc^ 


< 


‘^loc,k,g  —  V’fe— l,g— 1  (74) 

'^loc,k,q  —  0  OJJ-  Tfr} 


V'Jt.g/At  +  div{v'tpk,q)  -  t/Ai^k,q 

=  0  in 

'^k,q 

=  0  on  Too 

d‘^k,q 

9'^locyk^q 

dn 

~  dn 

(75) 


If  At  is  sufficiently  small,  we  can  apply  the  analysis  of  the  previous 
section  to  this  algorithm  and  we  conclude  that  V’fc.g  and  ^/oc,fc,g  converge 
linearly  to  zero.  Hence  the  sequences  and  are  Cauchy  sequences, 
which  converge  linearly  to  the  unique  solutions  and  of  the  implicit 
scheme.  This  guarantees  the  convergence  of  the  above  fixed  point  algorithm. 
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Appendix 


The  main  result  of  this  section  relies  on  the  notion  of  a  contact  set.  If  « 
is  a  continuous  arbitrary  function  on  fi,  the  upper  contact  set,  denoted  r+ 
or  r+,  is  the  sub-set  of  fi,  defined  by 


T'^  =  {y  e  3p(y)  €  lR”such  that  u(x)  <  u(y)  +  p  •  (x  -  y)yx  e 

(76) 

We  see  that  u  is  a  concave  function  on  fl  iff  r+  =  fi.  When  u  €  (7^(12) 
we  must  have  p  =  Du{y)  in  the  relation  (76).  In  addition,  when  u  € 
the  Hessian  matrix  D^u  =  [Diju]  is  negative  on  r+.  In  general,  r+  is  closed 
in  SI. 

If  M  is  a  continuous  arbitrary  function  on  SI,  we  define  the  "normal 
mapping”  xiv)  =  Xu{y)  at  point  y  e  SI  hy 

X(y)  =  {pe  IR”,  u{x)  <  u(y)  +  p.{x  -y)'ix  e  fl}.  (77) 

We  can  see  that  x{y)  is  non  empty  iff  y  €  r+.  In  addition  when  u  €  C^(Sl), 
we  have  x(2/)  =  J^u(y)  on  r+;  in  other  words  x  is  the  gradient  field  of  u  on 

r+. 

As  a  particular  case  of  the  Bahelman- Alexandrov  ([8]  and  [9])  maximum 
principle,  we  have  under  the  above  notation. 

Lemma  .1  For  u  €  C'^{Sl)  n  C°{Cl),  we  have  : 

sup u  <  supdQU  + 

“  nwn 

with  d  the  diameter  of  SI  and  Wn  the  volume  of  a  unit  sphere  in  R”. 

For  further  details  see  [12]. 

We  now  proceed  to  the  proof  of  Theorem  6.1,  by  following  the  steps  of 
[12].  We  take  B  =  Hi(0)  and  the  general  caae  will  be  deduced  by  considering 
the  coordinate  transform,  x  x  =  (x  —  y)  /2R. 

We  will  begin,  in  first  step,  by  showing  this  result  for  u  €  C^{Sl)  n 
W2’"(f2)  and  then  in  a  second  step  we  will  deduce  the  result  for  u  e  W2’”(fi). 
Step  1; 
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We  suppose  that  u  €  n  For  >  1,  we  consider  the  cut 

off  function  77  defined  by 

v(i)  =  (1  -  Ixff. 


By  differentiation,  we  obtain 


At?  =  -2pXi(l  -  \xff  \ 


Dijri  =  -2P6ij{l  -  \xff  ^  +  4p(p  -  l)xiXj{l  -  \x\'^f  ^ 


By  setting 


V  =  7]U, 


we  then  obtain 


a^^DijV  =  rja^^bijU  +  2a‘^  DirjDjU  +  ua^^Dijrj 

>  rj{f  —  VbiU  —  cu)  +  2a*^  biTfbjU  +  ua^^bijf]. 


Let  F'*'  =  r+  be  the  upper  contact  set  v,  in  the  sphere  B  ;  we  have  : 

u  >  0  on  F"^. 

If  X  €  dB  such  that  p.{x  —  y)  <  Owe  indeed  have  v{x)  =  0.  Consequently 

+  p-i^  -y)>  ‘o(x)  =  0- 

Moreover,  using  the  concavity  of  u  on  F"*",  we  can  estimate  the  following 
quantity  : 

\I)u\  =  (l/T?)|X)t;  —  UjC>t?|. 


Indeed, 


<  (i/>?)(| 


bv\  +  u 


bri\) 


<  (1 

<  2(1  +  y5)7?-i/^u. 


In  that  way,  we  have  on  F"*"  the  following  inequality  : 
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-a^^bijv  <  {(I6/32  +  27;/?)A77-2/^+ 

2/3|%“^/^  +  c}v  +  77I/I. 

Since  c  <  0,  we  deduce  the  inequality 

-a^^bijv  <  {(16/3^  +  277/?)A77"2/^  +  (78) 

2/3|6|7/"^/^}t;  +  77I/I 


<  C1A77  ‘^^^v  +  \f\, 

with  Cl  =  c(n,  /?,  7,  S)  independent  of  c. 

Consequently,  by  applying  Lemma  6.1  on  B,  we  obtain,  for  /?  >  2  : 


nWn 


By  using  the  relation  (2),  it  comes 

^  (— +  (— T7ir)(T)ll/lln,B 

B  nwr{  nw^  A 

<  + (i/A)ii/ii„g) 

<  c.d(||,-"*’<>-"ll.,4  +  (1/A)||/||,3) 

<  c,(i((s»p»+)>-2/'>||(ti+)W||„^g  +  (1/A)||/||,  g), 

where  ci  is  a  constant  depending  only  on  n,/3,7  and  6.  Here,  d  is  the 
diameter  of  B{d  =  2). 

By  using  the  Young  inequality  under  the  form 


ab  <  ea’’  + 


for  g  =  (1  —  2/j3)  ^  and  r  =  /3/2,  we  have 

<  esupv'^  +  Ve  >  0. 

By  taking  s  = - j  and  plugging  in  our  inequality  on  u,  we  obtain  : 

2cia 
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supv  <  (1/2)  sup v"'"  +  (1/2)^  (79) 

B 

+(cid7A)ll/iL,B- 

We  want  to  prove  the  theorem  for  all  p  >  0.  We  wiE  treat  separately 
the  cases  p  <n  and  p>  n. 

If  p  <  n,  we  set  ^  =  2n/p.  In  this  case  we  have 


Plugging  this  in  our  inequality  on  v,  we  obtain  : 

(l/2)sup„  <  (l/2)>-»=(c,i)«/"||(.+)||,3  +  (ci<i/A)||/||„p. 

B 

Consequently,  we  obtain  the  following  inequality  ; 

sap  t.  <  C2{(  I  (u+yylo  +  (d72A)||/||„,p}. 

B 

On  the  sphere  Bi/2(0),  the  cut  off  function  satisfies 

i/n  <  (1/2/- 

It  follows,  then 

sup  u  <  sup  {v/rj) 

Bi/2(0)  *61/2(0) 

<  2^supt;. 

B 

Finally  we  end  up  at  the  desired  estimate 

sap  u  <  C3{(  L  (n+rf/”  +  (<i/2A)||/||  p}. 

Bi/2(o)  Jb 

for  u  in  W^’”(fl)  (1  C^ip).  The  constant  C3  above  depend  only  on  n,  /3, 7  and 
6,  but  is  independent  of  c. 
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On  tlie  other  hand  li  p>  we  have  : 


2n/l5  <  p,  V/3  >  2. 

Then,  it  follows  (by  assuming  /?  >  2) 

But 


and  therefore,  by  processing  as  before,  we  obtain  the  desired  estimate 
sup  «  <  C4{(  f  +  (<i/2A)||/||  j} 

for  u^xn  n  C^(f2).  The  constant  C4  above  depends  only  on  n,/3,7 

and  Sj  but  is  independent  of  c. 

Transformation  x  x. 


By  construction,  Dij  =  thus  A  =  and  6  =  SB?,  In 

addition,  we  have  |B|  =  Wn{2RY‘  and 

Written  in  term  of  x,  the  last  inequality  becomes 

iZ  “  -  W 

with  C4  a  function  of  n,j,6  =  6R^  and  p.  This  is  the  desired  estimate  for 

u  6  n 

step  2: 

Now,  let  u  €  By  density,  let  (um)  be  a  sequence  of  functions 

of  C^(B),  converging  towards  u  in  The  injection  of  B)  in 

C°(B)  is  continuous,  consequently  (um)  converges  uniformly  towards  u  in 
B.  We  have 


LUm  =  L(u.m  —  u)  +  Lu 

>  f  +  L{Um  —  u). 
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By  setting,  fm  =  L{um  -  u),  we  observe  by  construction  that  fm  con¬ 
verges  towards  0  in  As  Um  €  D  C^(il)  and  fm  —  f  fm 

in  tbe  estimate  (3)  is  valid  also  for  Um,  so  that  we  have 

sup  U,u  <  C«e{(l4  /  «)")''”  +  7  I|/|U,b}-  (80) 

Bjiiy)  \^\  Jb  ^ 

Using  previous  results  and  taking  the  limit,  we  have  : 

sup  U<Cte{(y^  /  (ti+r)'/"+f|i/|kB}. 

SR(y)  l■®l  ^ 

m 

Observe  also  that  by  replacing  u  by  -u,  the  theorem  can  be  extended 
easily  to  the  case  of  super  solutions  and  solutions  of  the  equation  : 

Lu  =  /. 
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